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Abstract. The 1+3 covariant equations, embedded in an extended tetrad formalism 
and describing a space-time with an arbitrary energy-momentum distribution, are 
reconsidered. It is shown that, provided the 1+3 sphtting is performed with respect 
to a generic timehke congruence with tangent vector u, the Einstein field equations 
can be regarded as the integrability conditions for the Jacobi and Bianchi equations 
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together with the Ricci equations for u. The same conclusion holds for a generic null 
congruence in the Newman-Penrose framework. 
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. 1. Introduction 

CN 
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I Although the Einstein field equations 

Rij - ^Rgij = Tij (1) 

have a deceptively simple appearance, when written out in terms of coordinates and the 
components of the metric, they form an unwieldy set of second order non-linear partial 
differential equations, which is one of the reasons why, during the past half century 
and in several areas of gravity research (explicit construction and classification of exact 
solutions [T71 HB], cosmological perturbations and gravitational waves ([HI HI [TTl I5U]). 
numerical relativity ([ISl |13l [Ml EH 129]), fundamental aspects (O [3l |28l El 12]), focus 
has been shifting towards tetrad formulations of the theory. In this paper I consider the 
so called covariant 1+3 formalism and the related orthonormal tetrad formalism, which 
continue to play an important role in the context of cosmology, as well as the Newman- 
Penrose formalism. The structure of the governing equations of relativistic cosmology, 
including their consistency, in both formalisms has by now been discussed in a large 
number of papers[8l ESI UHl [20l [32l |35l |36] (for a recent review see the book by Ellis, 
Maartens and MacCallum [12]), while integrability of the equations in general tetrad 
formalisms has been discussed in detail by Edgar [6j and MacCallum|23|, following earlier 
work of Papapetrou|24l l25l 126] . There still remain some questionmarks regarding the 
redundancies present in these formalisms: while it is obvious that there are a number 



On the relation between the Einstein field equations and the Jacobi-Ricci-Bianchi system.2 



of algebraic interdependencies, the issue of differential relations (more particularly the 
possibility of obtaining certain equations as integrability conditions of others) is less 
clear. This paper aims to have a closer look at these differential interdependencies 
among the sets of Jacobi, Ricci and Bianchi equations on the one hand and the field 
equations on the other. As the commutator relations form a natural ingredient of such 
an investigation, the Jacobi equations will play a prominent role: I first repeat some 
basic facts about tetrad formalisms in section 2, closely following thereby paper [23] 
and pointing out some subtleties with regard to the Jacobi equations. In section 3 I 
show that the field equations turn out to be the integrability conditions of the Jacobi- 
Ricci-Bianchi set, provided the timelike congruence, with respect to which the 1+3 
splitting is carried out, has a non- vanishing acceleration which is not an eigenvector of 
aab + ojab- In section 4 I show that the same conclusion holds for a null-congruence k, 
provided that in a Newman-Penrose tetrad (m, m, £, k) the null-congruence k satisfies 
-K = ka-bk^m" ^ 0. 



2. Tetrad formalisms 



As usual space-time is represented by a pseudo-Riemannian manifold and all calculations 
apply to an open set of space-time. Notations and conventions are as in [18]. A 1+3 
splitting is accomplished locally by choosing a congruence of timelike worldlines with 
unit tangent vector field u. This congruence may be interpreted as the 4-velocity 
field of a family of observers, but no further restrictions apply: the congruence not 
necessarily consists of trajectories orthogonal to some family of spacelike hypersurfaces, 
neither needs the unit tangent vector field to be parallelly transported or be invariantly 
defined (although in most applications one of these situations might occur). A set of 
(smooth) basis vector fields is then constructed by erecting at each point a spatial triad 
orthogonal to (greek indices taking the values 1, 2, 3). Denoting the dual basis of 
one-forms as uj"" and the Levi-Civita connection as V, the connection coefficients and 
connection one- forms are given by V e^^b = '^'^ba^c and V^b = T^bc^^ respectively. The 
Cartan structure equations read then 

du"" = -T\ A oj^ (2) 

and 

dr\ + A r% = r\ , (3) 

R"-b being the curvature 2-forms. The integrability conditions for (l2]l3l) are given by the 
first and second Bianchi identities, d^a;*^ = and d^F^b = 0, or 

R\Au:'' = , (4) 
dR\ - R\ A T^b + T\ AR''b = . (5) 

Henceforth it will be assumed, as is usually done in tetrad formulations of general 
relativity, that the basis is rigid, in the sense that the metric components Qab are 
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constants, 

dgab = d{ea.eb) = , (6) 

with, for an orthonormal tetrad, gab = diag(— 1, +1, +1, +1). Raising and lowering 
tetrad indices with gab and its inverse and defining Tab = gac^^b & imphes Fab = ^lab] 
or 

^{ab)c = , (7) 

such that ([2]) completely defines the connection. 

The commutation coefficients 'y'^ab defined by 

[ea,eb]=Yabec , (8) 

or 

Yab = 2u^je[b\a] , (9) 

where e^* and a;°j are the coordinate component^ of the basis vectors and dual one- 
forms. The first Cartan equations ([2]) express that the connection is torsion-free and 
hence relate the connection and commutation coefficients by 

= 2r^[H , (10) 

which, for a rigid basis, by ([7]) is equivalent with 

Tcafe = -^{ibca + lacb " Icab) ■ (H) 

Introducing the components R^bcd of the curvature two-forms by 

i^^ = |i?^cd^^Aa;^ (12) 
the second Cartan equation becomes 

^ bed — >- 6d,c — J- 6c,d + J- ftd-L ec — J- bc>- ed — 1 cd^ be- U<jj 

The first Bianchi identity (jlj), which can also be written as 

Ra[bcd] = , (14) 

is then seen to be equivalent with the Jacobi identity for a triple of basis vectors e^, 

[e[a, [ei„ec]]] = , (15) 

or 

e[al\c] - l\\al\c\ = , (16) 

while the second Bianchi identity ([5]) can be written a^ 

R"'b[cd;e] = . (17) 

In a (pseudo-)Riemannian space the Riemann tensor also must satisfy 

Rabcd = —Rbacd = —Rabdc = Rcdab ■ (18) 

J a comma denoting a partial derivative or a directional derivative, depending on the index used 
§ a semi-colon denoting the covariant derivative with respect to V 
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There is some confusion about whether or not these symmetry conditions have to be 
imposed as extra conditions[23j: if the connection is defined in terms of the commutator 
coefficients by ( fTTl) and if the Riemann tensor is defined by ( fT3|) . then Rabcd = QaeR'^bcd is 
automatically anti-symmetric in the first and second pair of indices. Provided that the 
Jacobi equations ( !T6l) hold, this implies also Rabcd = Rcdab- This is what happens in a 
so-called mimmal tetrad formulatzon% 12^, where one regards the tetrad components 
Cq*, the 7^^^ and the matter fields as functions to be solved for, with the governing 
equations being the component forms and ( IT6|) of respectively the first Cartan and 
Jacobi equations and the field equations ([1]) together with the matter field equations. 
Note the crucial role of the Jacobi equations, which, via the symmetry conditions (JTSll . 
also imply that the Ricci tensor, defined by contraction of (fT3|) . 

Rbd = R'^bcd = ^e5^e,„-R'"d , (19) 

is symmetric. In addition the Einstein tensor is then divergence- free, such that the 
matter field equations must be compatible with T°'''.f, = 0. In this approach the 
remaining trace-free part of the curvature is given by the Weyl-tensor, which is defined 

by 

Cabcd = Rabcd — {ga[cRd]b + 9b[dRc]a) + lRga[c9d]b ■ (20) 

As noticed in [6] one can remove the tetrad components Cq* from the system and view 
the operators as the third set of variables alongside the ^'^^f^ and the matter fields. 
The governing equations are then the commutator relations ([8]), formally describing the 
behaviour of the 'variables', the Jacobi equations f|T6|) and the field equations ([1]) 
together with the matter field equations. 

Because of the special role played by the Weyl tensor in the construction of exact 
solutions or in the classification of space-time geometries, it has become customary to 
include the components of the Weyl tensor as extra variables in the previous system. 
In the case of an orthonormal tetrad formalism one decomposes then the Weyl tensor 
with respect to a timelike congruence defined hj u = Bq into its electric and magnetic 
components]^ 

C,,^'^ = 4{u[au^' + h^J')E,^ '1 + 2EabcU^'H''^^ + 2£^<^%[,//,]e , (21) 



or 



Eab = CabcdU^'u'' , (22) 
Hab = ^riac'^C.fMU'u" . (23) 

Provided the Jacobi equations hold and the Riemann tensor is defined by f[T^ . both 
Eab and Hab are then symmetric and trace-free. 



II The versions in the cited works differ shghtly. 

^ The Levi-Civita tensor is normalized such that (tetrad indices!) 770123 = —1; we also define 
Sabc = Tjabcd.u'^, whilc hab IS the projcctor in the instantaneous 3-spaces orthogonal to u, defined by 

hab = gab + UaUb- 
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Alternatively one can set up an extended tetrad formalism^l, in which the 
fundamental variables are 7^^^^ and the symmetric tensors Eab, Hat, Tab {Eab and Hab 
being also trace-free) and the matter fields. It is furthermore customary to split Tab as 

Tab = pUaUb + phab + 2q(^aUb) + T^ab , (24) 

where qa and Tiab = '^{ab) are orthogonal to Ua and tt^ = 0. In this extended formalism the 
Riemann tensor is denned by ( !20l) . using fin[2T]) to express it in terms of the variables Eab, 
Hab, T^ab, Qa, P and jQ. The symmetry and trace properties of these variables guarantee 
the conditions (HM ITSj) and the governing equations become the second Bianchi equations 
( IT7|) (written in terms of Eab, Hab, T^ab, Qa,P, P and their covariant derivatives) and (fT3|) . 
implying 1^. Note that ([13]) is now a set of partial differential equations for the 7*^60 
and not a definition for the Weyl tensor and that its contraction automatically yields the 
Einstein field equations ([1]). An equivalent approach is to implement (fT3l) by requiring 
that the Ricci equation holds, when applied to all four of the basis vector fields: 

«^^c;d - «'^d;c = R\dcW' , (25) 

again with appropriate substitutions for the Riemann tensor in terms of Eab, Hab and 

Tab- 

A third approach to implement (fT3l) (which is the one usually followed in 
cosmological applications, see for example §6.6 of [i2\) is to impose only the Ricci 
equation applied to u = eo, 

«";c;d - W'-Ac = R\dc (26) 

and to add the remaining field equations ([T]), together with the Bianchi equations (ITTl) 
and the Jacobi equations f lT6|) (which now not automatically hold, as only part of the 
Ricci equations is used). The variables of the resulting system (which obviously contains 
some redundancies) can be written exphcitly as oja, 0, aap, da, f^ap, ^a, P, P, la, 
Hap, Eai3, Hajj, with the kinematical quantities uj^ = ^SaP'y^^'^ , and being 
defined in the usual way by splitting Ua-b as 

Ua;b = -UaUb + CTab + \Ohab + OJab , (27) 
with aab = O-(ab), Cr^a = 0, Uab = UJlab], CTabU^ = ^abvf' = UaU"" = 0. 

= \r]°'^^'^Ubec ■ e.d is the local angular velocity, in the rest-frame of an observer with 
four-velocity m, of the triad Cq with respect to a set of Fermi-propagated axes (using 
MacCallum's convention |22]) and 71,^/3, eta are the Kundt-Schiicking-Behr variables[9] 
parametrizing the purely spatial commutation coefficients 7^/3^ ■ The commutation 
coefficients can be read off from 

[eo, e,] = li^eo - [\ebi + + ^^{u^ + fi'')) , (28) 

[e^, 6^3] = - 2eafi^u'e.fi + (^2a[„<5^] + e^fiin^^'^ (29) 

+ For example: i?i2i4 = -^31 - ^92, -R1212 = \p - ^'''33 - ... 
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and the connection one-forms read according! 

Tio = iiiuj^ + 9iuj^ + (ai2 + uj3)uj^ + ((T13 - ^2)^^ , (30) 
ri2 = -^^3^° + (nis - 02)0;^ + (n23 + ai)a;^ + ^(77,33 - - 7222)^1^^ -(31) 

In cosmological applications it is furthermore customary to define a fully projected 
covariant derivative V, having the property that for an arbitrary tensor 5" 

V,5^-"e.../ = K'h^ . . . h\hj . . . hf'SP-\...s,b . (32) 

This allows one -at least formally- to remove the Ua/s and Oq, variables from a subset 
of the equations ( nil6|ll7f26p . resulting in a set consisting of 4 Jacobi equations, 4 (Oa) 
Einstein field equations, 10 Ricci equations ( 126|) and the 20 Bianchi equations ( IT7|) and 
which is usually referred to as the "1+3 covariant equations". This system is quite el- 
egant, bears a strong analogy with the Maxwell equations |21| and forms the basis of a 
number of significant contributions to cosmology. It is an incomplete system and the 
fact that extra spatial information, representing the 12 remaining Jacobi equations and 
6 Einstein field equations, must be added in some way or another, has been reported 
on several occasions|10l [23l l35] . This incompleteness problem, in combination with the 
redundancies present in the fully extended system, was one of the motivations for the 
present investigation. It turns out that (see section [3]), when the timelike congruence 
corresponding to u is generic, adding the 6 spatial field equations is not necessary at 
all. In fact we will do more and show that generically the full set of field equations can 
be interpreted as integrability conditions for the system governed by ([51 [T6l) and (l26l) . 
This system (wherein the commutator relations ([8]) will always be implicitly assumed to 
be valid) will be referred to as the Jacobi-Ricci-Bianchi (JRB) system. The emphasis is 
on "generically", as it is clear that the mentioned result cannot possibly hold for every 
congruence in an arbitrary space-time: assuming for example a perfect fluid model with 
o'ab = = Ua = Eab = Hab = and taking u to be the fluid's timelike eigenvector]^ the 
only surviving JRB equations are 1) a single Bianchi equation (the conservation law for 
the matter density) 2) a single Ricci equation (the Raychaudhuri equation) and the 12 
Jacobi equations for ti^/s and a^, from which it is clearly impossible to derive the field 
equations. 

In the case of a null-tetrad formalism (Newman- Penrose), the situation is less 
complicated, as the equations are usually already presented in their 'extended' form, 
as sets of 'NP-equations' and 'Bianchi-equations', where the former is just the set of 
second Cartan equations (fT3l) . In section 4 I wil show how also the 'NP-set' can be split 
in Jacobi equations and Ricci equations for one of the null congruences (fe), with the 
remaining field equations being integrability conditions of the Jacobi-Ricci-Bianchi set. 

* For readability multiplets of equations, which can be obtained from each other by cychc permutation 
of the spatial indices, will henceforth be represented by a single equation. 

which, using Einstein's field equations, would then necessarily be a FLRW model 
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3. The JRB-system and the Einstein field equations in a 1+3 setting 

In order to have no misunderstanding about the meaning of the 'JRB-system' I first 
write out the equations explicitly. Where it is convenient, the expansion tensor 
©afe = o'ab + ^dhab is uscd in stead of the shear tensor and 9a = Qaa- Directional 
derivatives will be written as 9^ = e^. A hybrid notation will be used, where boldface 
symbols refer to objects with greek indices, while the 3-d operators ■ and x have their 
usual meaning, also when acting on non-tensorial objects: for example a ■ uj = OaOj'^, 
{d ■ n)a = dpn^ a, ■ ■ ■ (see also |33]). 

Jacobi equations 

Writing out the 16 Jacobi equations f lT6|) . one obtains 12 evolution equations. 



OqU) = — X u -\- ■ u + X u -\- cr ■ (jj + u; X fl — ^9uj , (33) 
doa =ld-fT-^dx{uj + fl)- \de - |(e + e{a + ii) 

+ {a-\u)x{ijj + n)-a-{a- \u) , (34) 
doni2 = \ (9i((T3i - ^2 - ^2) - ^2(0-23 + wi + J^i) + d-i^Qi - 6*1) 

+ 2(^11 + n22)o-i2 + 2(n22 - ^11) (^3 + ^3) - '^^12^^ 

+ Ml(0"31 - W2 - ^2) - M2(o-23 - Wl - - Usidl - ^2) 

+2n3i(cr23 + wi + ^^i) + 2n23(cr3i -a;2 -1^2)) , (35) 
donu = - di{uji + Qi) - 820^1 + 530-12 + (a + 2u) uj + d ■ Q, + ■ u 

- Ul{Ui + VLi) - U20-31 + U3Cri2 + nii{Oi - 6*2 - 6*3) 

+ 2^12 (cri2 + t^3 + ^3) + 2n3i(cr3i - U2 - O.2) (36) 

and 4 spatial 'divergence' equations, 

d ■ ijj = {2a + u) ■ u) , (37) 
dn= - dxa + 2@-uj + 2n-a + 2u}xn. (38) 



Ricci equations applied to u 

Writing out (p6l) and replacing, as explained above, the Riemann tensor by its 
decomposition in terms of Tab-, Eab and Hab, one obtains 18 independent equations, 3 of 
which are just the evolution equations fl33|) for the vorticity. The remaining set splits 
in 6 evolution equations for the expansion tensor, the trace being the Raychaudhuri 
equation, 

doO = _ 1 ^_ _ a^^a"^ + 2uaCo'^ + d ■ ii + ii ■ (ii - 2a) (39) 
and 

doOi = + diili -\e^ -\{p + 3p) - (Til ((Til + l^*) - (Tu - 4l + ^2 + ^3 

+ ul + U2{nsi - 02) - Uz{ni2 + 03) 
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+ 2(ai2^^3 - ^31^2) + Ittu - E^, , (40) 
doai2 = + d(iU2) - 0-12(^1 + ^2) - cr3ia23 + U1U2 - OJ1OJ2 + a(iM2) 

- \{nziui - n23U2) + ^usinii - 7222) 

+ (T3l^]l - CT23^^2 + {O2 - Oi)n3 + " ^12 ■ (41) 

There remain 9 equations for the spatial derivatives of the vorticity components, which 
are equivalent with the 9 'div cr' and 'curl cr' equations of the 1+3 covariant formalism, 

diU2 = dia-ii - 8361 + (o-23 - 2uiU2 + uJi){n3i - 02) + |a;3('^22 + "33 - '^ii) + 1^3 + ^12 
+ ^(^uinu + 3^22 - ^33) + (nu + a3)(^i - ^3) + 20-31(^23 - cti) , (42) 

diUs = 8261 - diai2 + (cr23 - uji){ni2 + 03) - 2uiUJ3 - 1^2(7233 + 7122 - Un) - \q2 + i/31 
+ |o-3i(72ii + 3^33 - 7122) + (7731 - a2){9i - O2) + 20-12(7723 + ai) , (43) 

diui = 930-12 - 52(^31 - ui^i + ^2(^2 + a2- 7731) + W3(m3 + a3 + 7712) + Hu 

+ 0-31 (r73i + 02) + CXi2(77i2 - ^3) " 27723^23 

+ i77n(3^i -6)- i(7722 - n33)(^2 - ^3) • (44) 

Note that with (jHj) also the Jacobi equation f l37|) becomes identically satisfied. 

Bianchi equations 

Writing out the Bianchi equations 0141) . one first obtains the 4 contracted equations, 
determining the evolution of p and qa, 

dop = -{p + p)e-d-q-2q-{u-a)- vr^^a"^ , (45) 
doHi = - dip - (p + p)ui - {d ■ Tz + Tz ■ (3a - u))i 

- qi{9 + 6*1) - 92(0-12 + a;3 - Vt^) - q^ia^i - a;2 + (^2) 

+ 7r3i77i2 - Vri2773i - (7r22 - 7r33)7723 + (^^22 " '^33)7r23 • (46) 

Then follow the sets of 10 'dot E', 'dot H' and 6 'div E', 'div H' equations, the appearance 
of which below is somewhat more complicated than the familiar (perfect fluid) one by 
the presence of the q^ and tt^^ terms: 

9o^l2 = - \d2qi - \ 5oVri2 + <92i^23 - <93^22 -\{uJ?, + (^12) {p + P) 

+ I {n22 - - 77n) 53 + | (-?t-23 - ai - ^i) ^2 - |m2?i 

+ {-Q2 + 2 0-31 - 2 U2) E23 + (0-23 + r^i - wi) ^31 + {-Q3 - 0-12 - ^3) Ell 

+ Q3E22 - 2 (ai2 + 0J3) E33 -{62 + 2 63) E12 

+ (772 - 2 a2 - 2 7731) H23 + (r7i2 - 03) -f^33 - (oi + + ^23) H31 + H11U3 

+ (a3 - 7712 - U3) H22 - \ (37711 + 7722 - ^33) ^12 

- |(7ril(0-i2 + U;3 + ^3) + ^27r23 + ^27ri2 - ^^37^"22 " (^1 " O-23 + ^\) T^3\) , (47) 
ao^ll = - l^oTTii - i^igi + 1(5393 + ^292) + 52^^31 - ^3/^12 + | (^ " S^i) (p + p) 

- I (2771 + ai) qi + l ((22 + 2u2 - Susi) 92 + | (03 + 377i2 + 2773) 93 
+ (cri2 + a;3 + 2 ^3) E12 - 2 £^23(723 + {(^31 - ^2 - 2 ^2) E31 
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-{02 + 2 63) + {02 - 63) S33 - iHunu + \ {7122 - n^z) {H22 - ^33) 

+ 2 i?23'^23 - (^^31 - 2 M2 + ^2) ^^31 - (^^12 + 2 M3 - 03) H12 

- \ (o'3i + 3a;2 + 6^2) 7r3i - \ {ai2 - - ^^z) t^vi + |7r23a-23 

- \ {26, + ^^3) VTll + ^{62- 63) 1122 , (48) 
doHi2 = d3{E22 + |P + ^TTn) - iaiTTgi - ^2^23 + | (0-3I - ^(^2) Ql " |?3^1 

+ i (3nii + 71,22 - ^^33) E12 + (ai + 7123 + ui) E31 + (2 7i3i -U2 + 2 02) E23 

+ 2 ^22^12 + (2 03 - U3) E33 + (-2 U3 + ni2 + 03) Ell 

+ (2 a3i -2uj2- ^2) H23 + ((T23 + - wi) //31 - (^2 + 2 ^3) H12 

+ (i:f22 - ^^^11)^3 + (0"12 + ^3) (-f^22 - i^33) " | + ^^12) (tTh - 7r33) 

+ i (^33 - nil - 37222) 7ri2 + (ai - 7123) t^z\ + |(a2 - 7231) 7^23 , (49) 
d^Hii = d3{Ei2 - jTTu) - d2{E3i - \ti3i) - qiui ^\{uj2- cr3i) q2 + \ (0-12 + ^^3) ^3 
+ {ni2 -03 + 2 U3) E12 - 2 E23n23 + (7131 + 02 - 2 U2) E31 + EuUu 

+ \ (^33 - nil - n22) E22 + \ ("22 - nil - ^^33) E33 - Hii{02 + 63) + H22O3 + i^33<92 
+ {ai2 +UJ3 + 2 9.3) H12 - 2 i^230-23 + (cr31 - ^2 - 2 9.2) H31 + \ (^22 - ^33) 7122 

- \ (3nii + 7733 - ^^22) TTii + 7r23ri23 " \ (^12 - ^3) 7ri2 - \ (tzsi + 02) 7r3i (50) 

and 

{d ■ E)3 = \d3P -\{d-TT)3 + luJiq2 - luJ2qi - \eq3 + i(e ■ 0)3 

+ 7131^23 - ^23^31 + ('^ll " ^22) -£^12 + (^22 " ^ll)'^12 + 3(-E ■ 0)3 

+ 0"l2(-f^22 — ^11) + (^1 — ^2) H12 + 0-31 i/23 — Cr23-f^31 " 3(i? ■ 0^)3 

+ |?^12(7r"22 - TTii) + i (77,11 - ^22) 7ri2 - |n237r3i + |r73i7r23 + f (tT ■ 0)3 , (51) 

{d ■ H)3 = \d2q1 - \d1q2 + uj3{p + + \a1q2 - \a2q1 + \{n- 0)3 

+ (6*2 — dl) E12 — cr3ii?23 + Cr23-E'31 + Cri2(-Eii — -E'22) + 3(£^ ■ Uj)3 

+ ^31^^23 - ^723^31 + (^1 " ?T-22) H12 + ni2{H22 - Hu) + 3(i/ • 0)3 

- 1(6*1 - O2) 7ri2 + |(T237r3i - |o-3i7r23 + 1^12 (vTn - 7r22) - \{-K ■ 0^)3 . (52) 

To sum up: the JRB-system consists of the following independent sets: 15 
Jacobi equations (12 evolution equations (133] - [36|) and 3 'div n' equations (138|) ). 15 
Ricci equations (6 evolution equations (1391 - I4T1) and 9 equations ( l42] - [44l) for the spatial 
derivatives of the vorticity and 20 Bianchi equations (14 evolution equations fl45|46|47f - 
I50|) and 6 'div E\ 'div H' equations fl5T]l52|) ). 

Einstein field equations 

We now write out the Einstein field equations, simplifying the right hand side of ([1]) 
by using all the algebraic information in the Jacobi and Ricci equations. The (Oa) 
components become then identically satisfied, while the (00) equation reduces to the 
trace of the {aa) field equations. We write the remaining 6 equations as yap = 0, with 

3^11 = 52(^^31 + 02) - 53(^^12 - %) + ^2^*3 - c^23^ + wi(t^i - 2^^l) - a^a" 
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+ |nii(2n22 + 2n33 - 3nu) + \{n22 - n^^f + 7123^ 
- ni2(ni2 - 203) - n3i(n3i + 202) + En - |p + Ittu 



(53) 



and 



3^12 = \di{n^i + ^2) - 1^2(^23 - oi) - \d'i{nii - n22) + iOi^2 + i^2^i - ^2^1 

- o"3i(T23 + ^'30-12 + ni2{nn + n22 - ^33) + 2 71317123 

- naiOi + 7i23a2 + 03(7111 - 7^22) - ^12 - |7ri2 • 

Integrability conditions for the JRB-system 



(54) 



The JRB-system provides us with expressions for all derivatives of the vorticity, in terms 
of spatial derivatives of the expansion tensor and the acceleration and, as is relatively 
easy to verify, the integrability conditions for these equations are identically satisfied 
under the Einstein field equations. That the field equations can also be viewed as 
integrabily conditions for the JRB-system, provided the u congruence is generic, can be 
seen by evaluating the following set of commutators, 

[do, di]{asi + 002) + [do, ds]ei + [ds, d,]ui , (55) 
[^0, 9i](cri2 + 0^3) - [do, d2]0i + [di, d2]ui , (56) 
[^1, 92](cr23 + uji) + [d2, d3]{ai2 - UJ3) + [ds, di]e2. (57) 

From the resulting expressions (and from similar ones obtained by cyclic permutation 
of the indices) all second order derivatives can be eliminated, leading to a homogeneous 
system of first order equations with the following simple structure: 

3^12 " 

3^23 = , (58) 

3^31 



iis 
-U3 Ui 

—Ml U2 



-U2 





diag(tii,ti2,M3) 



diag(7i2,«3,Mi; 



3^33 




3^11 








" 3^33 " 




3^11 




3^22 









7t3 
Ml 
M2 

-7i3 



-U2 
-(723 + Wi 





3^12 




3^23 




3^31 








3^12 




ill 




3^23 


= 







3^31 





(59) 



(60) 



Cr23 + 

-CT12 + W3 Cri2 + W3 

-(T31 + U2 0-31 + ^2 

— (731 — 0^2 02 — 03 Cri2 — <^3 

61 — 02 (731 — LJ2 — C"23 ~ ^1 

(T23 — cui — cri2— a;3 03 — 0i 

Aligning the 63 axis with the acceleration { iii = U2 - 

vanishing acceleration, = 3^23 = 3^3i = 3^ii = 3^22 



3^33 

3^11 

3^22 



+ 





3^12 




3^23 




3^31 



. (61) 

0), it is clear that, for non- 
= 0, which when substituted 
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in (IM]) implies (0-31 — ^2)3^33 = (a"23 + Wi)3^33 = 0. It follows that 3^33 = unless 
(7a3 + Wa3 = 0, meaning that u is an eigenvector of aab + oOab'- 

Theorem (1). If ii ^ is not an eigenvector of aab + ^ab, then the Einstein field 
equations are the integrability conditions of the JRB-system. 

Notice that this theorem is also an immediate consequence of the property that the 
curvature tensor for a given metric connection V is the unique tensor obeying 

• the symmetries (fTS!) 

• the first Bianchi equations (HM 

• the second Bianchi equations ( IT7|l 

• the Ricci equations (l26l) for some generic u congruence. 

A simple covariant proof of this uniqueness property is obtained by introducing the 
difference A of two such tensors and contracting ( ITTl) with m^, implying 

Aab[cdu''-e] = . (62) 

By (|26|) one has AabcdU'^ = and hence a further contraction of f l62|) with u'^ implies 
Aabcd.u'^ = 0. Defining then a tetrad with bq = u and 63 = -ii, it follows that the only 
possible non-zero component of A is given by j4i212 (= ^2121 = ~^i22i = "^2112)- 
Choosing in ( 162|) (acde) = (1123) or (2213) gives then 

(63) 



0, i.e. unless ii is an eigenvector of 



^2121^^^3 = A212M^;3 = 0, 

which indeed implies A = 0, unless u\^c.Ua]-bU^ 

(^ab + OJab- 

A similar property — not involving the acceleration of the time-like congruence — 
can be obtained by considering the integrability conditions of the second Bianchi 
equations: eliminating the second order derivatives from the commutator relations 
[do, d,]Hn + [do, d2]Hu + [do, d3]H,s-l[d2, ds]p-[d,, di]E^2-[di, d2]E,i-[d2, d3]E^, 
and [^0, di]Eu + [<9o, a2]^i2 + [^o, <93]^i3 - l[do, d3]p + [ds, dijH^s + [di, 82^^^ + 
[92: d^lH^i again results in a homogeneous system for the y^/^ variables: 



,(64) 



~2^23 





2^23 




3^33 




" -X31 


X22 — X33 


X12 




3^12 


X31 


-X31 







3^11 




2^23 


-X12 


X33 — Xii 




3^23 





X12 


-X12 




3^22 




2^11 ~ 2^22 


X31 


— X23 




3^31 





1^23 

with Xn 



-J21 



= E„ 



J12 


—J'i2 

^■Kab and Jc 





3^33 






3^11 






3^22 





'-ab — J^ab — 2'^ab anu Jab — Hab 

Generically the corresponding 6x6 determinant is different from zero and the system 
will again have only the zero-solution. Except for a few particular cases it is not easy 
however to provide a simple geometric characterisation of the non-vanishing of this 
determinant. It is clear that the left hand side of (IMI) is zero in an X eigenframe, such 



J22 — J\\ 

— Jz2 
J3I 

IVabcdQ^W^- 



-Jl3 
J\2 

Jzi, ~ J22 



J2Z 

Jll- 

-J2I 



J, 



33 





3^12 




3^23 




3^31 



(65) 
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that 3^12 = 3^23 = 3^31 = when X is non-degenerate. If I and H commute, a common 
eigen- frame can therefore be constructed in which clearly also = 3^22 = 3^33 = 
provided qiq2qz 7^ 0. Hence we obtain 

Theorem (2). // [X, i3"] = 0, X^p is not degenerate and if q is not parallel to one of 
the eigenblades of H, then the Einstein field equations are the integrability conditions 
of the JRB-system. 

On the other hand it is easy to verify that the Einstein field equations cannot be 
obtained from ( !M|65|) when for example q = 0, as then the matrix in the left hand side 
of ( l65l) is singular. 

4. The case of a null congruence 

When the congruence, with respect to which the Ricci equations are constructed, is null, 
one can construct a Newman-Penrose tetrad (m, m, i, k) with k = u: the equations 
are then usually presented |27[ [T8] in their 'already extended' form, as the set of 'NP- 
equations' — these being just the second Cartan equations (fT3l) — and the second Bianchi 
equations. We will refer to the two sets by their numbering in [18j as N Pi-N Pig and 
Bi-Bii, with the latter being just the set of integrability conditions of the former. One 
may again isolate from the NP-equations the subsets of Jacobi and Ricci equations and 
ask whether the remaining part of the NP-equations (the 'remaining' field equations) 
is obtainable as integrability conditions of the Jacobi-Ricci-Bianchi set. An explicit 
evaluation of f lT5|) shows that the Jacobi-part of the NP-equations is expressed by the 
16 independent equations 

Im - NP^, m\i - NPii, WK2 + iVT^ - NP^2 - NP^, 
NPe - NPe + NPs - NI\, NP^q - WPj, NP^^ + NPs, 
2iVF4 - ^ - NP^, 2NP, - NPs + NPu, 

2NP,, -Jm- NPTs, 2NP,s + NPg - NP^s . (66) 
The Ricci equations f l25|) for k on the other hand give rise to 18 independent equations 
NP,, NP2, NPs, NPn, iVPie, NP^r, NP^ + NP,, 

iVPis - NPis,lm + NP^,NPr2-NPu , (67) 
which together with fl66|) reduce the remaining field equations to a set of 6 equations, 

iVPio, NPis, NPu + NPr2, NPu + NPr4 (68) 

(or any set equivalent with it under (16611671) ). Considering the linear combinations of 
commutators [S, 6]k,+ [6, D]p-[6, D]a, [A, D]p+[6, A]k-[S,D]t, [A, D]a-[6, D]j + 
[S, A]e and [A, D]f3 — [6, D]^ + [6, A]e one obtains then, using (I66|67p . the equations 
kNPis, kNPu, kNPu - kNPio and 

D{NPi3) + - (p - 2e)A^Pi3 - 2(7r + t)NPi2 + 2kNPio , 

from which, provided n ^ 0, (168 p readily follows: 
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Theorem (3). If n = —ka-hm"'k^ ^ the Newman- Penrose equations NPiq,NPi3, 
NP12 + NP12 and NP14 + NP14 are the integrability conditions of the remaining Newman- 
Penrose equations and Bianchi equations. 

5. Conclusion 

Embedding tlie 1+3 covariant equations in an extended tetrad formalism [12] leads to 
redundancies, forcing one to make a choice among the sets of Jacobi, Ricci, Bianchi 
equations and Einstein field equations (whereby with 'Ricci equations' we mean the full 
set of Ricci equations applied to the tangent vectorfield of a single time-like congruence). 
It is shown that a minimal set of equations can consist of 15 Jacobi and 15 Ricci 
equations, together with the 20 Bianchi equations, in which the Riemann tensor is 
defined in terms of trace-free and symmetric tensors Eab-, Hab and a symmetric tensor 
Tab- The Einstein field equations arise then as the integrability conditions for this set, 
if one chooses the time-like congruence such that the acceleration is not zero and is 
no eigenvector of aab + ^ab) a condition which obviously is violated for pressure-free 
matter. It remains to be seen whether this alternative view of the Einstein equations, 
as integrability conditions of an extended system, has any useful applications at all, for 
example in the area of numerical relativity. In the case of a null congruence n = fc an 
analoguous result is obtained provided that the associated Newman- Penrose coefficient 
K = — fcajfem'^fc* is non- vanishing. 
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